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One of the most striking astrophysical phenomena today is the existence of cosmic ray particles with 
energies in excess of 10^® eV. While their presence has been confirmed by a number of experiments, 
it is not clear where and how these particles are accelerated to these energies and how they travel 
astronomical distances without substantial energy loss. We are entering an exciting new era in cosmic 
ray physics, with instruments now producing data of unprecedented quality and quantity to tackle the 
many open questions. This paper reviews the current experimental status of cosmic ray physics and 
summarizes recent results on the energy spectrum and arrival directions of ultra-high-energy cosmic 
rays. 


1 Introduction 


Cosmic ray particles were discovered almost 
one hundred years ago, and yet very little 
is known about the origin of the most ener¬ 
getic particles above and around 10^® eV, tra¬ 
ditionally referred to as “ultra-high-energy 
cosmic rays.” The measured spectrum of cos¬ 
mic rays extends beyond 10^° eV, 11 orders of 
magnitude greater than the equivalent rest 
mass of the proton. While the presence of 
particles at these energies has been confirmed 
by a number of experiments, it is not clear 
where and how these particles are accelerated 
to these energies and how they travel astro¬ 
nomical distances without substantial energy 
loss. Some indication comes from the en¬ 
ergy spectrum itself, which roughly follows 
a power law ® and is therefore “non- 
thermal.” The power law behavior and the 
“universal” spectral index indicate that the 
underlying acceleration mechanism could be 
Fermi acceleration^ at shock fronts, e.g. in 
Supernova Remnants (SNRs) and in the jets 
of Active Galactic Nuclei (AGNs). Regard¬ 
less of the actual acceleration process, it is 
clear that thermal emission processes cannot 
generate such energies. 

Several quantities accessible to experi¬ 
ment can help to reveal the sources of ultra- 
high-energy cosmic rays: the flux of cosmic 
rays; their chemical composition] and their 


arrival direction. Gharged cosmic ray pri¬ 
maries are subject to deflection in Galac¬ 
tic and intergalactic magnetic fields and do 
not necessarily point back to their sources. 
The strength and orientation of these fields 
is poorly known and estimates vary^’®, but 
their impact should decrease at the highest 
energies, above several times 10®® eV; here, 
cosmic ray astronomy might be possible. 

Since the cosmic ray flux drops almost 
three orders of magnitude for each energy 
decade, the flux at the highest energies is 
very low, about one particle per km^ per 
year above 5 x 10®® eV. Low statistics have 
historically plagued the field; the total pub¬ 
lished number of events above 4.0 x 10®® eV 
is still less than 100, and drawing conclu¬ 
sions on the basis of such a small data set 
has proved rather perilous. However, the 
experimental situation is finally improving; 
new instruments are now collecting data of 
unprecedented quality and quantity. Since 
1999, the High Resolution Fly’s Eye (HiRes) 
air fluorescence stereo detector in Utah has 
accumulated data with excellent angular res¬ 
olution. The HiRes data set has been used 
extensively over the last two years to search 
for small-scale anisotropies in the arrival di¬ 
rections of cosmic rays, correlations of cosmic 
rays with known astrophysical sources, and to 
study the composition of the primary cosmic 
ray flux. 
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Figure 1. Scheme of the different cosmic ray detec¬ 
tion techniques (Auger Collaboration). 


In the southern hemisphere, the world’s 
largest detector for cosmic radiation, the 
Pierre Auger Observatory in Argentina, is 
nearing completion, and first results on the 
energy spectrum and the arrival direction dis¬ 
tribution of cosmic rays have recently been 
published. 

In this paper, I will review some re¬ 
cent developments in ultra-high-energy cos¬ 
mic ray physics. The paper is organized as 
follows: Section 2 gives a short review of the 
experimental techniques and the current ma¬ 
jor instruments in the field. I will then dis¬ 
cuss new results from the Auger and HiRes 
experiments on the energy spectrum (Section 
3) and the arrival direction distribution (Sec¬ 
tion 4). Concluding remarks follow in Section 
5. 


2 Experimental Techniques 

Because the flux at ultra-high energies is 
small, experiments need a large detector vol¬ 
ume. Consequently, detectors have to be 
earth-bound, and the primary cosmic ray 
particles can not be observed directly. Pri¬ 
maries interact in the upper atmosphere and 
induce extensive air showers with on the or¬ 
der of 10^° particles for a 10^® eV primary. 
The properties of the original cosmic ray par¬ 
ticle, such as arrival direction and energy, 
have to be inferred from the observed prop¬ 
erties of the extensive air shower. 

There are two different techniques to 
study cosmic ray air showers at ultra-high- 
energies. Both are shown schematically in 
Fig. 1. Ground arrays sample the shower 
front arriving on the Earth’s surface with an 
array of particle detectors, for example scin¬ 
tillation counters or water Cherenkov detec¬ 
tors. The arrival direction of the air shower 
and the primary cosmic ray particle is recon¬ 
structed from the differences in trigger times 
for individual detectors as the narrow shower 
front passes. The advantage of ground arrays 
is their near 100 % duty cycle and the robust¬ 
ness of the detectors. A disadvantage is that 
ground arrays sample the shower at one alti¬ 
tude only and do not record the development 
of the shower in the atmosphere. Moreover, 
the sampling density is typically very small. 

The classic example of a pure ground ar¬ 
ray is the AGASA (Akeno Giant Air Shower 
Array) experiment, which operated in Japan 
from 1984 to 2003. In its final stage, the ar¬ 
ray consisted of 111 scintillation counters of 
2.2 m^ area each on a 1 km spacing, leading 
to a total area of about 100 km^. 

Apart from the cascade of secondary par¬ 
ticles, air showers also produce Cherenkov 
light and fluorescence light. The latter is pro¬ 
duced when the particles in the air shower 
cascade excite air molecules, which fluoresce 
in the UV. Air fluorescence detectors mea¬ 
sure this light with photomultiplier cameras 
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Table 1. Comparison of Instruments. Comments: above 10^^ eV, 10% duty cycle, January 2004 — June 

2004, 3) error bars are strongly asymmetric 


Experiment 

Operation 

Aperture 
[km^ sr] 

Exposure 
[km^ sr yr] 

Angular Resolution 
(68%) 

AE/E 

[%] 

AGASA 

1984-2003 

~ 250 

1620 

2.5° 

25 

HiRes mono 

1997 - 

10000^^ 

5000 

> 2.5°’®^ 

25 

HiRes stereo 

2000 - 

10000^) 

3400 

O 

o 

15 

Auger (under 
construction) 

2004 - 

7400 

1750 (SD)2) 

2.0° - 0.9° (SD) 
0.6° (FD) 

10 


that observe the night sky. The shower is 
observed by a succession of tubes and re¬ 
constructed using the photomultiplier tim¬ 
ing and pulseheight information. Air fluo¬ 
rescence detectors can only operate on clear, 
moonless nights with good atmospheric con¬ 
ditions, so the duty cycle is about 10 %; how¬ 
ever, they observe the shower development in 
the atmosphere and provide us with a nearly 
calorimetric energy estimate. In addition, the 
instantaneous detector volume is rather large, 
of order 10000 km^ sr at 10^° eV. 

The High Resolution Fly’s Eye (HiRes) 
experiment in Utah is a stereo air fluores¬ 
cence detector with 2 sites roughly 13 km 
apart. Each site is made up of several tele¬ 
scope units monitoring different parts of the 
night sky. With 22 (42) telescopes at the 
first (second) site, the full detector covers 
about 360° (336°) in azimuth and 3° — 16.5° 
(3° — 30°) in elevation above horizon. Each 
telescope consists of a mirror with an area of 
about 5m^ for light collection and a cluster 
of 256 hexagonal photomultiplier tubes in the 
focal plane. 

The HiRes air fluorescence detector can 
operate in “monocular mode,” with air show¬ 
ers only observed from one site, or in “stereo¬ 
scopic mode,” with the same shower ob¬ 
served by both detectors simultaneously. The 
monocular operation suffers from poor angu¬ 
lar and energy resolution. With only one 
“eye,” the shower-detector-plane {i.e. the 
plane that contains the shower and the de¬ 
tector (see Fig. 1)) can be reconstructed with 


high accuracy. Unfortunately, the position of 
the shower within that plane, determined us¬ 
ing the photomultiplier trigger times, is am¬ 
biguous. Stereo viewing of the shower with 
two detectors breaks the ambiguity and leads 
to an excellent angular resolution of order 
0.5°. HiRes has been taking data in monocu¬ 
lar mode since 1996 and in stereo mode since 
December 1999. 

In summary, stereo air fluorescence de¬ 
tectors have excellent angular resolution and 
give a nearly calorimetric energy determina¬ 
tion, while ground arrays have the advan¬ 
tage of a relatively straightforward determi¬ 
nation of the detector aperture. Obviously, 
the best detector is a detector that com¬ 
bines the two techniques. The Pierre Auger 
Observatory, currently under construction in 
Malargue, Argentina, and scheduled to be 
completed in 2006, is such a hybrid detector, 
combining both a ground array and fluores¬ 
cence detectors. 

The Auger surface detector array^ will 
eventually comprise 1600 detector stations 
with 1500 m separation, spread over a total 
area of 3000 km^. One year of Auger data- 
taking will therefore correspond to about 
30 AGASA years. Each detector station 
is a light-tight 11000 liter tank filled with 
pure water. Three 9-inch photomultipliers in 
the tank measure the Cherenkov light from 
shower particles crossing the tank. The sta¬ 
tions are self-contained and work on solar 
power. The surface detector array is com¬ 
plemented by 4 fluorescence stations® with 6 
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telescopes each. The field of view of a sin¬ 
gle telescope covers 30° in azimuth and 28.6° 
in elevation, adding up to a total field of 
view of 180° X 28.6° for each site. Each tele¬ 
scope consists of a spherical mirror of area 
3.5 m X 3.5 m and a focal surface with 440 
hexagonal (45 mm diameter) photomultipli¬ 
ers. Because of the large field of view of each 
telescope, a Schmidt optics is used. 

While still under construction, the Auger 
experiment has recorded data with a growing 
detector since January 2004. The total ex¬ 
posure of the surface detector array already 
exceeds the total AGASA exposure, making 
Auger a competitive experiment. 

Table 1 gives an overview of experimen¬ 
tal parameters for AGASA, Auger, HiRes 
monocular and HiRes stereo, including the 
aperture of the instruments and the exposure 
used in publications of recent results. 


3 Results 

Until a few years ago, the world data set of 
ultra-high-energy cosmic rays was dominated 
by data recorded with the AGASA air shower 
array. In a number of important publica¬ 
tions, the AGASA group has described sev¬ 
eral exciting and controversial results, includ¬ 
ing the shape of the energy spectrum above 
10^® eV, studies of possible anisotropies in the 
arrival directions above 4.0 x 10^® eV, and 
an excess of cosmic ray flux from the Galac¬ 
tic center region around 10^® eV. Due to the 
small number of events, some of these results 
have a small statistical significance. For sev¬ 
eral years, statistically independent data sets 
to support or refute these findings were not 
available. Only recently, with the start of 
HiRes stereo data-taking in 1999 and the be¬ 
ginning of operations at the Pierre Auger Ob¬ 
servatory in January 2004, are we reaching 
the point where we can study these topics 
with new data. 


3.1 Energy Spectrum 

Proton primaries above 5.0 x 10^® eV inter¬ 
act with the 2.7 K microwave background via 
photo-pion production, losing energy in each 
interaction until they eventually fall below 
the energy threshold. This so-called GZK 
effect, postulated by Greisen® and indepen¬ 
dently by Zatsepin and Kuzmin^ shortly after 
the discovery of the microwave background, is 
expected to lead to a rapid fall-off of the cos¬ 
mic ray energy spectrum above this energy, 
but it has not been experimentally confirmed 
at this point. 

The AGASA group claimed®’® that the 
GZK suppression is not observed, raising 
questions as to the nature of the primary 
particles or even the particle physics in¬ 
volved. More recently, the HiRes monocular 
spectrum^®’^^ has been interpreted as being 
in agreement with a GZK suppression. While 
the disagreement between the AGASA and 
the HiRes mono result has received a lot of 
attention, it has also been pointed out that 
the two spectra actually agree reasonable well 
if systematic errors are taken seriously^®’^®. 

Fig. 2 depicts the situation. It shows the 
differential energy spectrum as reported by 
AGASA® and HiRes^®’^^ in monocular mode 
(the HiRes collaboration has not yet pub¬ 
lished a stereo spectrum). While only sta¬ 
tistical errors are shown, systematic errors of 
order 25 % (roughly what is reported by the 
experiments) would manifest themselves sim¬ 
ply by a shift of one bin in this plot. Keeping 
this in mind, both spectra agree quite well be¬ 
low 10^®'® eV. The disagreement at the high¬ 
est energies has been estimated to be of order 
2.0(7 12. 

Reversing the argument, it is actually 
quite surprising how well the spectra agree, 
given that ground arrays and fluorescence de¬ 
tectors have entirely different systematic er¬ 
rors in their estimate of the shower energy. 
AGASA has given a detailed account of their 
energy resolution®. Ground arrays determine 
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the energy of a shower from the measured 
particle density S at some (fixed) distance 
from the shower core^®. The optimal distance 
is the one where the spread in particle den¬ 
sity is minimal, and this optimal distance de¬ 
pends on the detector geometry (600 m for 
AGASA® and 1000m for Auger^®). The cor¬ 
relation between the particle density and the 
energy is then established using simulated 
data. This of course leads to a strong mass 
and model dependence. 

Air fluorescence detectors image the 
shower development in the atmosphere and 
obtain a nearly calorimetric signal. However, 
the Earth’s atmosphere is a tricky calorime¬ 
ter. It is not only highly inhomogeneous, 
but also changes on short time scales. The 
atmosphere needs to be monitored continu¬ 
ously in order to be able to correctly account 
for Rayleigh and Mie scattering. In addition, 
there are systematic errors from the fluores¬ 
cence yield, the absolute calibration, and the 
accounting for “missing energy,” i.e. energy 
that goes into particles that do not contribute 
to air fluorescence. 

The hybrid nature of the Pierre Auger 
Observatory will enable us to understand 
whether there are intrinsic difficulties with 
one of the two methods that contribute to 
the disagreement in the two spectra. At 
this time, the Auger experiment has not 
yet accumulated the necessary statistics to 
make a definitive statement on the shape 
of the energy spectrum at the highest ener¬ 
gies. The collaboration has, however, already 
published a first estimate of the spectrum. 
Fig. 3 compares the Auger spectrum to both 
the AGASA and HiRes 1 monocular measure¬ 
ment. The spectrum is based on 3525 events 
with energies above 3 x 10^® eV data taken 
with the ground array between January 2004 
and June 2005. At these energies, the ground 
array is fully efficient and the total exposure 
can be inferred from the total running time 
(where the growing size of the detector is ac¬ 
counted for) and the geometrical aperture of 


the ground array. The energy of the show¬ 
ers is established using the subset of hybrid 
events in the data sample, i.e. the subset of 
events that are also detected in one of the flu¬ 
orescence detectors. These events are used to 
derive a relationship between shower energy 
(as determined by the fluorescence detector) 
and the ground parameter S (as measured by 
the ground array) 

By using the statistics and the exposure 
of the ground array and the nearly calorimet¬ 
ric energy estimate of the fluorescence de¬ 
tectors, it is ensured that the Auger spec¬ 
trum does not rely strongly on either sim¬ 
ulations or assumptions about the chemical 
composition of the cosmic ray flux. Unlike 
the AGASA and HiRes spectra, it is there¬ 
fore nearly model-independent. 

As shown in Fig. 3, the Auger spectrum 
agrees quite well with the HiRes measure¬ 
ment at this point. Gonsidering the large 
systematic uncertainties, there is also little 
disagreement with AGASA. For this first es¬ 
timate of the spectrum, systematic errors 
range from about 30 % at 3 x 10^® eV to 50 % 
at 10^® eV, of which about 25% stem from 
total systematic uncertainties in the fluores¬ 
cence detector energy measurements. Errors 
will soon become considerably smaller with 
larger statistics and the application of more 
sophisticated reconstruction methods. 

Two additional aspects of the Auger 
spectrum should be stressed. 

(1) While the spectrum does not show 
any events above 10^® eV, the Auger ex¬ 
periment has detected an event above that 
energy^^. The event does not enter the en¬ 
ergy spectrum because its core falls outside 
the surface detector array, but the event is 
well-reconstructed and passes all other qual¬ 
ity cuts. This means that all major cosmic 
ray experiments have now confirmed the ex¬ 
istence of particles above 10^® eV. 

(2) As described above, the energy scale 
of the spectrum is normalized to the fluores¬ 
cence detector. When simulations are used 
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Figure 2. Ultra-high-energy cosmic ray flux as a function of energy as measured by the AG AS A experiment^ 
and the HiRes detectors in monocular mode, HiResl^^ and HiRes2^^. 
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Figure 3. Intensity of the ultra—high-energy cosmic ray flux as a function of energy measured by the 
AGASA experiment^ and HiRes 1 in monocular mode^®, compared to a first estimate by the Pierre Auger 
Observatory^^. 
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to determine the shower energy from sur¬ 
face detector data alone, the energies are sys¬ 
tematically higher by at least 25%. Ener¬ 
gies determined in this way depend on the 
shower simulation code, the hadronic model, 
and the assumed composition of the cosmic 
ray flux. The difference shows that many 
systematics need to be addressed. Auger is 
in a unique position to study these system¬ 
atic uncertainties and will measure the spec¬ 
trum in the southern hemisphere accurately 
in coming years. There is no reason to be¬ 
lieve that the energy spectra in the northern 
sky and the southern sky are identical, al¬ 
though differences are probably subtle and 
will require a larger instrument in the north¬ 
ern hemisphere. 

3.2 Search for Small-Scale Clustering 
and Point Sources 

A direct way to search for the sources of 
ultra-high energy cosmic rays is to study 
their arrival direction distribution. Astron¬ 
omy with charged particles of course faces 
a serious problem. At low energies. Galac¬ 
tic and intergalactic magnetic fields will bend 
the particle’s trajectories sufficiently to ren¬ 
der the direction information useless. How¬ 
ever, the Larmor radius increases with en¬ 
ergy, and cosmic ray astronomy may be pos¬ 
sible above some energy threshold. But little 
is known about intergalactic magnetic fields, 
so it is not straightforward to determine an 
optimal energy cut for arrival direction stud¬ 
ies. Choosing an energy threshold too low 
means that deflections destroy any correla¬ 
tion, but too high a threshold also weakens 
the statistical power of the data set. Given 
these uncertainties, an a priori optimal choice 
for the energy threshold and the angular sep¬ 
aration for clustering searches does not exist. 

AGASA claimed a significant amount of 
clustering in the arrival direction distribution 
as early as 1996^® and has updated these re¬ 
sults frequently^®’^®’^^’^^. The current claim 


is that the 5 “doublets” and 1 “triplet” ob¬ 
served in the set of 57 events with energies 
above 4.0 x 10^® eV have a probability of less 
than 0.1 % to occur by chance in an AGASA 
data set of this size. This is an extremely im¬ 
portant result, but its validity has been ques¬ 
tioned on statistical grounds. The problems 
stem from the fact that the data set used to 
evaluate the chance probability includes the 
data used for formulating the clustering hy¬ 
pothesis in the first place; parameters like 
the energy threshold and the angular scale 
that defines a “cluster” of arrival directions 
were not defined a priori. The probability of 
0.1% has therefore little meaning, and anal¬ 
yses that test the hypothesis only with the 
part of the data set that was recorded after 
the hypothesis was formulated find a much 
higher chance probability of around 8 % , 

making the result insignificant. 

The AGASA clustering claim has been 
tested with the statistically independent 
HiRes stereo data set. The HiRes col¬ 
laboration has published several papers de¬ 
scribing searches for deviation from isotropy, 
including the calculation of the standard 
two-point correlation function^'^, a two-point 
correlation scan^®, and a search for point 
sources using an unbinned maximum likeli¬ 
hood ratio test^®’^"^. None of these searches 
have uncovered any deviation from isotropy, 
and no evidence for statistically significant 
point sources in the HiRes or the combined 
HiRes/AGASA data sets above 4.0 x 10^® eV 
or 1.0 X 10^® eV were found^®’^^ (although 
a point source has been claimed by other 
authors^®’^®). With an overall exposure that 
exceeds the AGASA exposure and an angular 
resolution that is 4 to 10 times sharper than 
AGASA’s, HiRes stereo finds no small-scale 
anisotropy at any energy threshold above 
10^® eV and any angular scale out to 5°. 
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3.3 Correlation with BL Lac Objects 

There has recently been marginal evidence 
that a small fraction of the ultra-high-energy 
cosmic ray flux originates from BL Lacer- 
tae Objects^°’^^. BL Lacs are a subclass of 
blazars, which are active galaxies in which 
the jet axis happens to point almost directly 
along the line of sight. The EGRET instru¬ 
ment on board the Compton Gamma Ray 
Observatory (CGRO) has firmly established 
blazars as sources of high energy 7 -rays above 
100MeV^^, and several BL Lac objects have 
been observed at TeV energies with ground- 
based air Cherenkov telescopes^^. High en¬ 
ergy 7 -rays could be by-products of electro¬ 
magnetic cascades from energy losses asso¬ 
ciated with the acceleration of ultra-high- 
energy cosmic rays and their propagation in 
intergalactic space^^’^®. 

The history of claims for a correlation 
between BL Lacs and ultra-high-energy cos¬ 
mic rays is rather convoluted. Correlation 
claims based on data recorded with AGASA 
and the Yakutsk array have been published 
since 2001 ^®’^”^’^®, but a problematic aspect 
of the claims is the procedure used to es¬ 
tablish correlations and evaluate their sta¬ 
tistical significance. The authors explicitly 
tuned their selection criteria to assemble cat¬ 
alogs showing a maximum correlation with 
arrival directions of cosmic rays above some 
energy. The statistical significance quoted 
is therefore meaningless, and claims of BL 
Lac correlations have been criticized on these 
grounds®®’'*^. In some cases it has been shown 
that statistically independent data sets do 
not confirm the correlations^^. 

Recently, the HiRes collaboration pub¬ 
lished its own analysis of the data®^, us¬ 
ing an unbinned maximum likelihood ra¬ 
tio test which accounts for the individual 
point spread function of each event. The 
HiRes stereo data does not confirm any of 
the previous claims based on AGASA and 
Yakutsk data, in spite of its larger statisti¬ 


cal power. It does, however, verify a recent 
analysis®^ of correlations between published 
HiRes stereo events above 10^® eV and a sub¬ 
set of confirmed BL Lacs from the 10th Veron 
Gatalog^^. This subset®® contains 157 con¬ 
firmed BL Lacs from the catalog with optical 
magnitude m < 18. Since the cuts®® used to 
isolate this signal are not a priori^ the correla¬ 
tion needs to be confirmed with statistically 
independent data before any claims can be 
made. 

The correlations are strongest at energies 
around 10^® eV, where magnetic fields are ex¬ 
pected to sufficiently scramble the arrival di¬ 
rections of charged particles, yet the corre¬ 
lations are on the scale of the detector an¬ 
gular resolution. This would suggest neutral 
cosmic ray primaries for these events, or at 
least that the primaries were neutral during 
significant portions of their journey through 
Galactic and extragalactic magnetic fields. If 
verified with future data, the correlation with 
BL Lac objects would be the first evidence 
for an extragalactic origin of the highest en¬ 
ergy cosmic rays, and a first indication that 
at least a fraction of this flux originates in 
known astrophysical objects. 

3.4 Galactic Center 

One of the most interesting regions of the sky 
to search for an excess of ultra-high-energy 
cosmic rays is around the Galactic center. 
This region is a natural site for cosmic ray 
acceleration. It harbors a black hole of mass 
2.6 X 10® solar masses whose position is con¬ 
sistent with the radio source Sagittarius A*. 
Hour-scale X-ray and rapid IR flaring indi¬ 
cate the presence of an active nucleus with 
low bolometric luminosity. In addition, this 
crowded part of the sky contains a dense clus¬ 
ter of stars, stellar remnants, and the super¬ 
nova remnant Sgr A East. The Galactic cen¬ 
ter is now also established as a source of TeV 
7 -radiation"‘®’"‘"‘’"‘®. 

The AGASA experiment reported an ex- 
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Figure 4. Auger skymap in the vicinity of the Galactic center^®, using surface detector data. Shown are the 
coverage map (upper left) and the significance maps for the data smoothed using three different resolutions: 
the Auger resolution (upper right), the large integration radius used by the AGASA collaboration for their 
claim of a Galactic center source (lower right), and the resolution of the SUGAR experiment (lower left). The 
Galactic center is indicated by a cross, the regions of excess for AGASA and SUGAR are indicated by the red 
circles. The dashed line indicates the limit of the AGASA field of view. 


cess with a statistical significance of order 
4(7 near the Galactic center'^®’'^^. The ex¬ 
cess is observed in a narrow energy band from 
1017-9 eV to eV only. Due to its loca¬ 

tion in the northern hemisphere and a zenith 
angle cut of 60°, AGASA’s field of view cuts 
off roughly 5° north of the Galactic center, so 
the center itself, at right ascension a = 266.4° 
and declination S = —28.9°, is outside the 
field of view. The AGASA excess is found 
roughly around a = 280°, S = —16°, so it is 
offset considerably from the location of the 
Galactic center. To produce this excess, the 
event density is integrated over a circle with 
a radius of 20°. Several other “beam sizes” 
were tried, but this integration radius was 
found to maximize the signal. 

The chance probability in the AGASA 
publication is a posteriori. Since the anal¬ 
ysis has several “tunable” parameters like 
the energy bin and the integration radius, 
the true chance probability cannot be de¬ 


rived from the AGASA data set itself. The 
AGASA result was, however, supported by 
a re-analysis of archival data taken with the 
SUGAR array‘d®, a cosmic ray experiment 
that operated in Australia between 1968 and 
1979. Unlike AGASA, SUGAR operated 
from a location with good visibility of the 
Galactic center region. The SUGAR data 
showed an excess of about 2.9 tr in the energy 
bin 10^’^ ® eV to 10^® ® eV, roughly the same 
energy bin as AGASA, but at a = 274° and 
5 = —22°, so offset both from the location of 
the AGASA excess and the Galactic center. 
Furthermore, the SUGAR excess was consis¬ 
tent with a point source, indicating neutral 
primaries. 

In such a context, it is interesting to 
point out that neutron primaries are a viable 
hypothesis, as they could travel undeflected. 
The neutron hypothesis would also explain 
the narrow energy range of the signal. In 
an amazing coincidence, the neutron decay 
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length at 10^® eV roughly corresponds to the 
distance between us and the Galactic center, 
about 8.5 kpc. 

Like SUGAR, the Auger detector has 
good visibility of the Galactic center region. 
The Auger coverage map in the vicinity of 
the Galactic center is smooth and shows 
no strong variations (see upper left part of 
Fig. 4). The Auger group has searched for 
an excess in the Galactic center region in the 
data set taken from 1 January 2004 through 
5 June 2005^®, both in the surface detec¬ 
tor data (angular resolution 1.5°) and hy¬ 
brid data (0.6°). The event statistics for 
the surface detector data are already larger 
than those of the two previous experiments 
(3 times AGASA, 10 times SUGAR), but 
the AGASA excess occurs in an energy range 
where Auger is not fully efficient (> 30 % for 
protons, >50% for iron). However, even in 
the worst case (a proton signal on an iron- 
dominated background) a 5.2 cr excess is ex¬ 
pected if the AGASA source is real. 

Fig. 4 shows the results of an analysis 
using surface detector data taken between 
January 2004 and June 2005. The figure 
shows the map smoothed using different an¬ 
gular resolutions, corresponding to the Auger 
resolution (upper right), the AGASA sig¬ 
nal (lower right), and the SUGAR resolution 
(lower left). No significant excess is found 
in all cases. The Auger group also studied 
neighboring energy bins to account for a pos¬ 
sible systematic difference in the energy cal¬ 
ibration, but no excess is found in any sce¬ 
nario. For the Galactic center itself, a 95 % 
confidence level upper limit for the flux from 
a point source is derived^® which excludes the 
neutron source scenario suggested to explain 
the previous claims^^’®®. 

4 Concluding Remarks 

The most important and most fascinating re¬ 
sult at this point is that cosmic ray particles 
with energies above 10^® eV exist - their exis¬ 


tence has been confirmed by all experiments, 
regardless of the experimental technique, and 
regardless of whether the GZK suppression is 
present or absent in the data. 

With current statistics, the cosmic ray 
sky is remarkably isotropic. One possible rea¬ 
son for this could be magnetic smearing. This 
effect should be smallest at the highest ener¬ 
gies, above 10^® eV, but until now the num¬ 
ber of events has been too small to draw any 
conclusions. In another possible scenario, we 
could be dealing with many sources, each cur¬ 
rently contributing at most one or two events, 
and again, more data will help to eventually 
resolve the strongest ones. 

The Pierre Auger Observatory will dra¬ 
matically increase the high-energy sample 
size over the next several years. Auger is still 
under construction (to be completed in 2006), 
but the collaboration has already published 
first results, among them a first estimate of 
the cosmic ray energy spectrum. With its lo¬ 
cation in the southern hemisphere, its hybrid 
design and unprecedented size. Auger is in an 
excellent position to answer definitively the 
questions left open by the previous genera¬ 
tion of cosmic ray experiments, and to make 
the discoveries that will challenge the next 
one. 
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5 Discussion 

V. Gavasinni (Univ. di Pisa/INFN): If the 
cosmic ray excess from the Galactic center 
is due to neutrons, would you expect also a 
neutrino flux coming from the neutron decays 
which could be measured by large volume neu¬ 
trino detectors ? 

Yes. Neutrons with energy less than EeV will 
decay in flight and produce an antineutrino 
flux above TeV that could be detected by the 
next generation neutrino telescopes such as 
IceCube. The expected event rate per year 
above 1 TeV for IceGube was estimated to be 
about 20. See reference [51] for details. 

T. Greenshaw (Liverpool University): Does 
the simultaneous observation of the fluores¬ 
cence and surface signals at Auger give ad¬ 


ditional information on the composition of 
UHE cosmic rays and, if so, do you have any 
preliminary results ? 

There is a wealth of information in the com¬ 
bined measurement, but studies are ongoing 
and there are no published results from Auger 
yet. The HiRes collaboration has published 
results based on the stereo data. See refer¬ 
ence [52] for details. 

J. Lefrancois (LAL, Orsay): You showed a 
plot with the Auger and AGASA energy spec¬ 
tra. Gan the different rate of almost a factor 
of three be explained by energy calibration ? 
The rate differences between AGASA, HiRes 
mono, and Auger indeed suggest problems 
in the energy calibration. Further evidence 
comes from the first analysis of the Auger 
spectrum. Here, it has already been shown 
that there is a discrepancy in energies if the 
spectrum has its energy scale normalized to 
the fluorescence detector or if simulations are 
used to determine the shower energy from 
surface detector data alone. This might in¬ 
dicate problems with the shower simulation 
code, the hadronic model, and the assumed 
composition of the cosmic ray flux. The dif¬ 
ference shows that many systematics need to 
be addressed. 




